subroutine masscalc(x,mass)
! Subrotina para calcular a massa da partícula no sistema
! Pode ser modificada de acordo com o sistema
! Requer apenas a entrada da posição x e retorna o valor da
! massa naquele ponto do sistema

! *** Sem Variáveis Implícitas
	implicit none
! *** Declaração de Constantes e Variáveis
! * Constantes
	double precision, parameter :: nx = 10001
! * Variáveis - Argumentos Burros
	double precision, dimension(nx) :: x
	double precision, dimension(0:nx+1) :: mass
! * Variáveis - Argumentos Não-Burros
	integer i
	double precision dx, Ry, a0, am, am1, am2, ec, pi
	double complex :: img, dt

	common / constants / dx, dt, img, Ry, a0, am, am1, am2, ec, pi

do i = 1,nx
	mass(i) = am1
	if(x(i).GT.(0.D0)) mass(i) = am2
enddo

mass(0) = mass(1)
mass(nx+1) = mass(nx)

return
end subroutine masscalc
